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Abstract 

We analyze the role played by temperature in QCD3 by means of a dimensional inter- 
polating approach. Pure gauge QCD3 is defined on a strip of finite width L, which acts 
as an interpolating parameter between two and three dimensions. A two-dimensional 
effective theory can be constructed for small enough widths giving the same longitu- 
dinal physics as QCD3. Explicit calculations of T-dependent QCD3 observables can 
thus be performed. The generation of a deconfinig phase transition, absent in QCD2, 
is proven through an exact calculation of the electric or Debye mass at high T. Low 
and high T behaviors of relevant thermodynamic functions are also worked out. An 
accurate estimate of the critical temperature is given and its evolution with L is studied 
in detail. 
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1 Introduction 



In the non-perturbative nature of non-abelian gauge theories in three and four dimensions 
stems the complexity of the confinement phenomenon. On the contrary confinement shows 
up as an outgrowth of one dimensional space kinematics in 2D non-abelian gauge theories. 
This is not a striking fact if one considers the importance of kinematical constraints in 2D 
Quantum Field Theories, e.g. in conformal field theories. It is therefore far from obvious 
what kind of connection links these confinement mechanisms in different dimensions, not 
even if such connection exists. 

In a previous work |J we have studied pure gauge QCD 3 on a strip of width L. This 
theory is defined on a 3D (2+1) space-time manifold having pure gauge QCD 2 (L — > 0) and 
QCD3 (L — > 00) as extreme limits. An important property of this dimensional approach is 
that it allows us to go beyond the pure QCD 2 results in a controlled way (the e expansion, 
e being the dimensionless parameter e = g^L 1 ^ 2 ). No matter how small e is, one is dealing 
with a real 3D theory provided e 7^ 0. This fact liberates the theory from the strong 2D 
constraints we previously mentioned, and at the same time provides us with a panoply of 
new physical phenomena. The explicit solution of the problem is obtained by the analysis of 
a 2D effective theory equivalent to the original 3D one in the regime of very small widths: the 
theory of heavy adjoint matter coupled to 2D gluons in a gauge invariant way. This theory 
possesses a complete glueball spectrum which can be exactly computed for small values of 
the strip width. Interesting results of various versions of this model can be found in recent 
literature JJ. 

In the dimensional approach the singularity and simplicity of the theory's line limit 
(QCD 2 ) arises due the absolute dominance of pure longitudinal gluons (2D gluons) when 
the width vanishes. All transverse effects are excluded (brute force!) by the requirement 
of pure 2D kinematics. As soon as we allow the strip to acquire a width, transverse gluons 
start to play a role. Their appearance transforms the physical content of the theory in 
a highly non trivial way. They are responsible for a rich new glueball spectrum having 
relevant consequences in the thermodynamics of the theory. As it will be shown throughout 
this paper, the new physics associated to transverse gluons is behind the existence of a 



deconfining phase transition at high temperature. 

In order to understand the connection between the transverse dynamics and temperature, 
we will calculate the free energy of two external static charges in pure gauge QCD 2 - 

Let us begin by defining our fields in a 2D space-time with a compactified Euclidean time 
direction (as usual in the imaginary time formalism ||) 

a a (x,t) = a° a (x) + J2<(x)e l2nnt/P (1) 

We will work in the ai = gauge, where the gluon self couplings disappear. The charges are 
represented as external static color sources, which we will take in the fundamental represen- 
tation and forming a color singlet state (they behave like heavy quarks in a meson). 

Since the sources are static only the gauge field zero mode will couple to them. Thus we 
just need to perform a Gaussian integral to obtain the free energy of the pair: 
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(2) 



The external source corresponding to two static color charges in the singlet channel is 
given by j b = g 2 [Q b S(x - x Q ) + Q b 5(x - xq)] where -Q b Q b = C N = (N 2 - 1)/2N. 

Therefore we obtain a linear growth of the energy with R, the string tension being 
a = (gf|/2)Cjv- Besides there is no temperature dependence in the free energy. This is an 
interesting feature that is worth analyzing in the framework of the strip formalism. In a 
strip with no width, in which there are no transverse degrees of freedom, the stringlike con- 
figuration joining both charges is the only one contributing to the free energy. The stringlike 
configuration is nothing but the classical solution of static pure gauge QCD 2 . Longitudinal 
perturbations around this classical solution are forbidden because of the presence of the two 
infinitely heavy sources. Transverse perturbations are not allowed by the kinematical con- 
straint induced by the absence of a width. Therefore, as it is manifest in the free energy 
calculation (0), the system cannot jump to any other configuration at any temperature. The 



probability of the system to be in the classical stringlike configuration is always one. In other 
words the entropy of the system is zero for all temperatures and thus the free energy becomes 
T-independent S = ~(dF/dT)y = 0. Dynamics is frozen and the system lies always in a 
confining phase. 

The previous thermodynamic argument shows already why pure longitudinal gluon physics 
cannot "unfreeze" the 2D string joining the static quarks in the meson. It is clear now the 
relevance of transverse degrees of freedom in this picture. They are responsible of the ther- 
modynamics of the quark- ant iquark pair. The possibility of giving some small degree of 
transverse dynamics to the system in a controlled way becomes then crucial to understand 
the role of temperature in the previous scenario. This is precisely the purpose of our con- 
tribution. By using an effective 2D theory valid for narrow strips we will try to unveil the 
interesting connection between traversal dynamics, confinement and temperature in pure 
gauge QCD 3 . 

2 Characterization of confinement 

The v.e.v. of the Polyakov loop, defined as, 



is an order parameter of confinement depending on temperature. This is so because it is 
related to the free energy of an isolated static color pointlike charge through the expression 



<0|tr c £(0)|0> T 



(3) 



where 




(4) 
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e -(l/T)F Q (T Q ;T) =<0 | tr ^ (0) | 0>r (5) 

where Tq is the color isospin of the static charge (we refer to SU(2) -in general we would 
need the set of Casimirs characterizing the representation of SU(N)-), and the matrix g is 



taken in the representation labeled by Tq. If the v.e.v. of the Polyakov is zero we are in 
a confining regime. The energy to bring the isolated charge from infinity to x = would 
be infinite. On the contrary if it is not zero, we might liberate it (sending it to infinity) by 
supplying a finite amount of energy to the system: deconfinement occurs! The conventional 
wisdom about this order parameter relates it to the discrete symmetry associated to the 
center of the gauge group (Zn). Since under gauge transformations the Polyakov operator 
tr c g is color invariant up to an element of this discrete group tr c g— »Z7vtr c g, its v.e.v. is also 
an order parameter of the Zn symmetry. 

Nevertheless we are not going to adopt this approach to characterize confinement in our 
problem. Since we describe heavy sources with quark quantum numbers, the operator g(x) 
is a SU (N) matrix in the fundamental representation. As a SU (N) matrix, one can perform 
not only local vector transformation upon it, but more general ones. In particular, we can 
distinguish left and right color indices and rotate them independently with SU (N) matrices. 
In other words, we can perform SU(N) L ® SU(N) R transformations on g(x) in the following 
way, 



g(x) -> U L (x)g(x)U R (x) 
U L (x) = e lTa <^ 

U R (x) = e^^W (6) 

If we choose a a L = a a R then U L = U R = U and we recover the usual color vector trans- 
formations under which (tr c g) is invariant (up to the center). However under axial color 
rotations, defined through the conditions a a L = —a a R — (3 a (Ul — U R — Ua) , the Polyakov 
operator transforms non-trivially and its trace is no longer invariant (tr c g— >ti c {U\g)). Con- 
sequently, the v.e.v. of the Polyakov operator becomes also an order parameter of color axial 
symmetry breaking. Ordinary symmetry breaking analysis holds. If the vacuum preserves 
the color axial symmetry then <ti c g> T = and thus confinement occurs. If not, this v.e.v. 
will acquire a non-zero value and the system will be in a deconfining phase. 

We investigate the confinement phase structure of pure gauge QCD 3 on the strip through 
the analysis of its behavior under such symmetry. The fundamental field in the study of the 



Polyakov loop is the gauge field zero mode. We are concerned with its transformation rule 
under Ua{N). For 577(2) it takes a specially simple form. If tt(x) stands for the modulus of 
the gauge field zero mode and /3 for that of the axial angle, then the referred transformation 
reads, 

tt(x) ^ tt(x) + 2(3 (7) 

Color vector transformations keep the modulus unchanged. Axial ones shift it by a 
constant. Moreover if we perform a global vector rotation in such a way the gauge field zero 
mode at x = gets aligned along the 3rd direction of color isospin, the above transformation 
is equivalent to a shift in a 3 (0), (from now on a a stands for the zero mode field omitting the 
0- subscript) 

a 3 (0)^a 3 (0) + 2/? (8) 

This is possible because both the strip action and the operator ti c g are color singlets. 
Due to translation invariance, the v.e.v. of ti c g is independent of position. Thus Eq.@ 
would apply to any value of x. Notice that these transformations constitute the maximal 
abelian subgroup of the axial color non-abelian group Ua- 

We have found that we can characterize confinement by studying the behavior of the 
gauge field zero mode dynamics under global U% rotations, as defined by Eq.(^j). The vacuum 
properties of the a-field effective theory become crucial to grasp the confinement mechanism. 

The symmetry properties of the vacuum can be more easily understood in the Schrodinger 
formalism. To study the vacuum structure we need to know the behavior of the a field at 
low energies (long distances). The long distance action for the gauge field zero mode will be 
ruled by the terms with less number of derivatives, 



S[a a ] = Jdx -(1 + Z{a 2 )){d x a a f + V(a 2 ) + ... (9) 

where we choose V(0) = 0. Notice that Z and V depend only on the modulus of the gauge 
field due to color invariance. 



The previous Euclidean action is one dimensional. Thus we are really dealing with 
a Quantum Mechanical (QM) problem in the Euclidean "time" x. If we focus in the 
Schrodinger equation for the ground state of the QM counterpart of Eq.(^), we will find 
(a a ^q a , d x a a ^j^p a , p a = -ig|r) 

*o[g°] = E ^o[q a ] (10) 

The vacuum energy Eq is the minimum of all the eigenvalues of the QM Hamiltonian 
operator of the above equation. The Euclidean action (^) is positive defined, thus E > 0. 
Because of color invariance the wave functions of the QM problem depend only on the 
modulus of the gauge field zero mode ^/[g 2 ]. Therefore they are invariant under vector 
rotations, q — > UyqU v (recall the gauge field zero mode is no longer a gauge degree of 
freedom). We are free to choose a vector transformation Uy rotating q in such a way the 
color isospin vector lies always along the third component axis. That is, ty[q a ] = ty[(q 3 )u v ]. 
Let us see what happens if besides we demand the wave function to be invariant under U3 
transformations. This would be precisely the case of the vacuum wave function of the strip 
action in the confining phase, as we have seen above. A U3 symmetric wave function must 
satisfy, 

d 

— V[q 3 } = =>• * = const (11) 
oq 6 

Therefore the wave function of the U3 invariant vacuum is independent of the gauge field 
zero mode. If we plug the above constant vacuum wave function in the ground state equation 
Eq.fllOD, we find two new properties. First, a q dependent potential term is not compatible 
with such a solution. Thus V = (recall we have chosen V(0) = 0). Second, the ground 
state energy is zero. 

The Us invariant vacuum wave function (confining phase) is a constant whose energy 
eigenvalue is zero. Moreover the effective potential V cannot exist in this phase. The value 
of this constant defining the ground state wave function, Eq. (|TTD , is completely arbitrary 
and so the wave function is. Thus the ground state becomes degenerate. This degeneracy is 
nothing but that induced by the U3 symmetry. 



V 2 + V(q 2 ) 

2(l + Z(g 2 )) V " + {q ] 



Let us go into more detail at this point. A constant gauge field zero mode (a = const) 
minimizes the action (Q) where there is no potential V for symmetry requirements. Therefore 
each ground state wave function is generated by a constant field a configuration, and vice 
versa. We can move from one zero energy configuration to another one just by changing the 
value of the constant field. This is certainly what a U3 transformation does when acting on 
a constant a fielcf]. 

However if the theory develops a non-zero effective potential for the gauge field zero mode 
a, the previous property is no longer true. The presence of an effective potential prevents 
the vacuum wave function to be a constant. The U3 degeneracy is lost. In the broken phase 
the effective potential must be different from zero. Since spontaneous breaking of the U3 
symmetry cannot occur in such a low dimension (a dynamics is defined in D = 1) the only 
possibility for the vacuum to break it is through an explicit potential term (again this can 
be seen by studying the Schrodinger equation for the vacuum wave function Eq. (|10|) ) . 

Therefore one can use the effective potential for the gauge field zero mode itself as an 
order parameter of confinement. A zero effective potential implies we are in a confining 
phase. A non-zero effective potential indicates deconfinement occurs. 

We can check the previous pattern of confinement in pure gauge QCD 2 Eq.(||]). In this 
simple example the effective action for the ao is automatically obtained just by removing the 
non-zero mode fields a n . They do not couple either to the source nor to the zero mode, so they 
can be integrated out changing only the overall normalization constant. After integration 
the effective action for the gauge field zero mode contains only the kinetical term, which is 
certainly invariant under U% shifts, and no effective potential appears at any temperature. 
This agrees with the fact that pure gauge QCD 2 is always in a confining phase. 

3 The generation of the deconfining phase 

It has been proven |1| that a 2D effective action can be constructed to describe the physics 
of pure gauge QCD 3 on a very narrow strip. The effective action corresponds to adjoint 
1 Notice that for a constant field a U3 transformation is completely equivalent to a whole axial one XJ a- 
We can always choose the gauge field to be aligned along the 3rd direction for all x. 



scalar Chromodynamics, that is, to QCD2 coupled to color adjoint matter, 



^strip = n\fa, + + ^(e)0 2 ] (12) 

e being the dimensionless parameter e = g 2 L. The 2D coupling constant g 2 is related in 
turn to the 3D one, 5 , 2 = fi , |/e- In the regime where the above action is well-founded, e <C 1, 
adjoint matter is heavy since 

^ = f%(lA 2 ) = gln(lAr-°<x> (13) 

Eventually it decouples, leaving pure gauge QCD 2 as the effective 2D theory of the 
widthless strip. Notice that Eq.flT2|) entitles us to go beyond QCD 2 dynamics in a well 
defined way. The field is nothing but the zero mode of the transverse component of the 
gluon field ((f) ~ J L dc 2 A 2 ). It represents then a 3D transverse gluon moving on a very 
narrow strip. The new dynamics introduced by this transverse degree of freedom changes 
radically that of pure gauge QCD 2 . 3D transverse gluons (in the effective form of 2D heavy 
scalar adjoint matter) bind themselves in color singlets by means of the confining interaction 
supplied by the presence of longitudinal gluons. They make up the glueball spectrum of 
QCD 3 on a narrow strip |]. 

Apart of being interesting on its own, the existence of a non-trivial glueball spectrum has 
other consequences. At the very moment the system acquires new configurations to access 
to, thermodynamics ceases to be trivial. The partition function and all the thermodynamical 
properties derived from it become temperature dependent. 

In a 3D language, the action (^) reproduces the dynamics of transverse fluctuations 
around the classical, and unique, solution for e = (A^ = <5 M o«o + 0)- Unlike in pure gauge 
QCD 2 one may now excite kinematically the classical string configuration in the transverse 
direction by giving some energy to the system. Whether these thermodynamical fluctuations 
are able to break the longitudinal string confining the charges, or not, is a matter of the 
particular dynamics enclosed in the effective action (|T2]). 

In order to investigate the high temperature properties of QCD on the strip, we choose 
the a\ = gauge and select both the gauge and the scalar field to be periodic in the Euclidean 



time coordinate. We can expand them in terms of their ID Fourier modes (as in Eq.([T])) and 
obtain the following ID strip action at temperature T, 



l r. ri 



S(T) = f Jdx {-dl + f{e) + ul)5 n ^ n ,5 ab + (M^f" + \ (14) 



where, 



u n = 27rnT 

(M 2 ) n a r' = gh abC t bdf <»*-n-n'-n» (^) 

The (j) integration can be formally done yielding the effective action for the gauge field 
modes, 

S eff [a n ] = i Jdp ^Tr [\n(p 2 + /i 2 + u 2 + M x + M 2 )] (16) 

Now the trace runs over both color and the mode index n. The calculation of that trace 
is impossible in the most general case. However expression ([TED is useful to determine the 
limits in which a well defined approximation can be established. 

We are concerned with the possibility that the theory undergoes a deconfining transition 
at high temperature. According to what was discussed in the last section, we need to know 
whether the dynamics is able to generate a non-zero effective potential for ao or not. At 
zero T such a term is not allowed since a is a truly gauge degree of freedom in that case. 
Local gauge invariance prevents its existence. At finite T the gauge field zero mode becomes 
a covariant object under gauge transformations. As a consequence we cannot longer resort 
to gauge invariance to protect the theory from getting a non-zero effective potential. This 
possibility becomes therefore a dynamical issue. 

At very high T and for very narrow strips, for which adjoint matter is very heavy, the 
following inequalities hold, 



T>^(e) > g 2 



c^ + /i 2 ( e ) > /i 2 (e)>s 2 2 (17) 

The previous regime legitimates a perturbative expansion of the logarithm in the effective 
action in terms of the gauge field modes (|i"6l). The coupling constant dependence arises 
through the Mi and M 2 operators Eq. flT5|) . Their contributions are small when compared to 
those arising from the non-interaction terms. 

After performing the gauge field expansion keeping only terms up to order g%, one obtains 
a typical finite T 1-loop calculation for the mass of the /-mode (see Fig.l), 



<nUc,T) = dT\Y,l^ - 2 , - 2 , „, -£/•£ ,_, , ,., ""t"f' + " — 2 r ) (18) 



fer) = 9ir {?/l ^ + ^+^ -?/i 



The first, and very important, property of the above mass is that it is zero for all non-zero 
modes. Every non-zero mode is massless to this order of the approximation. We will point 
out the relevance of this feature in our conclusions. 

On the contrary, the gauge field zero mode acquires a non-zero mass which is computed 
systematically under the conditions given by the inequalities in Eq . (|17|) . Therefore at high 
temperatures the electric, or Debye, mass is given by 



m 2 el (e,T) 
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(19) 



T 'ix 

Moreover one can prove that, within the regime determined by the conditions (|i"7D, there 
are no corrections to the derivative dependent part of the effective action (fLED (they are 
suppressed as powers of g\j 'fx 2 ). In this way we get an exact high temperature result for the 
zero mode effective action of QCD on the strip, 



SefflM = -fj dx 



\al{-^al + \ml{e,T)alal 



92,4 



(20) 



Our question about whether transverse gluon fluctuations, even if small, would be able to 
break the tube flux structure at high enough temperature has now a clear answer. According 



to our discussion in the previous section the appearance of the mass term in Eq.(^) is enough 
to ensure the breaking of the U3 axial symmetry and, as a consequence, the generation of 
a non- vanishing confinement order parameter (<tr c g>^ 0). We conclude then that QCD 3 
defined on a strip experiences a deconfinig phase transition at high temperature. 

But we can go further, for we can explicitly check this result by calculating the Wilson 
loop (W(T,R)) in an analogous way as we did in Eq.(H). We only have to substitute the 
linear ID propagator of pure gauge QCD2 by the exponentially decreasing one derived from 
the action (p0|). The functional integration is immediate because is quadratic in oq, 



-lnW(T,R) = jpFqq(T, R) "g" 1 jdx \ft( x ) ( \ ^ j ,%(x) 

= I-f-p.^)*?*'' (21) 
T m e i m e iT 

In the a± = gauge the Wilson loop can be easily related to the Polyakov operator, 
W(T,R) —< tr c [g(R)g'(0)] >. Because of cluster decomposition, our above result certainly 
implies, 



W(T,R) ^ 2 <ti c g(0)> 2 => <tr c g>= e 2T ^ (22) 

In agreement with the general arguments given in the previous section. Notice how the 
absence of a mass is the responsible of the vanishing of the v.e.v. of the Polyakov operator. 
The strong IR behavior of the gauge field zero mode turns the ordinary analytical expansion 
of the exponential operator into an essential singularity. 



4 Low temperature glueballs 

The theory at zero temperature exhibits a confining behavior. The Wilson loop fulfills an 
area law at long distances in the same way as pure gauge QCD2- The only difference appears 
in the exact value of the string tension, which is now renormalized by quantum transverse 
fluctuations (0-loops) 0, ||. The spectrum of the theory at small e can be found in reference 
The low lying spectrum is calculated explicitly and is formed by heavy glueballs, made 



up of two heavy constituent transverse gluons. The radial wave functions are proportional 
to Airy functions with energies given by, 



( 



a 



2 



) 



1/3 



M r (e) 



2// fl (e) + e r 



(23) 



(— e r are the zeroes of the Airy function and a is the — string tension). The renor- 
malized scalar mass fin is obtained by using a self-consistent treatment of the IR divergences 



The low temperature behavior of the partition function will be ruled by the lowest mass 
states of the physical Fock space. They are eigenstates of the Hamiltonian operator at 
zero temperature. Because confinement occurs at zero temperature the only physical states 
contributing to the partition function are glueballs (singlet — boundstates). In the most 
general case (even when vacuum excitations and non-planarity effects are suppressed, e.g. in 
the iV— >oo) the eigenstates of the full Hamiltonian will not coincide with those obtained by 
solving the bound state equation of the — system Eq . (|23"D . This is so because glueballs 
can interact among them in an effective way [Q and a further mass renormalization can 



In our case we are going to consider glueball states as almost-free particles, exposed only 
to contact interactions. Therefore our Hamiltonian eigenstates will be given by the solutions 
of the two body — bound state equation. This approximation is completely justified 
for small enough e since for very heavy glueballs, as it is the case, the glueball interaction 
can be neglected. The OBEP (One Boson Exchange Potential) is extremely short range 
(~ 1/Mqb <S Rgb , M GB being the glueball mass and Rgb the glueball radius), as one can 
easily estimate using the strip results. The glueball radius can be computed using the known 
radial wave functions. If we call R = (yUR<7)~ 1//3 then the glueball radius can be obtained 
from the mean squared radius of the two particle system, 



0. 




(24) 
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=<x 2 > 1/2 - 1.707510 i? ( 
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(25) 



Therefore in the small e limit, the correlation length of the glueball exchange is negligible 
when compared to the typical radius of the glueball because 



n ^l/(M G M) e ^(^^0 (26) 

K GB 

On the other hand the typical radius of a glueball Rqb decreases logarithmically as the 
strip width tends to zero, R GB ~R = (/i^tx)- 1 / 3 ~ (ln(l/e 2 ))- 1/3 ^ 0. We conclude that 
glueballs are extended small objects in the low e regime. 

Consequently our system is made up of small quasi-free glueballs in the small width 
regime. Moreover since they are also very heavy, a non-relativistic approach (NR) is likewise 
appropriate at low temperature. For this reason the NR free glueball gas model becomes 
an excellent approximation to the real theory at low temperature. Since the dimensionless 
width e is a free parameter, we can play with it in order to make our approach as good as 
possible. In this way the narrower the strip is the better the free gas approximation is. 

The thermodynamics of the NR gas of glueballs is described by a partition function 
which factorizes in individual partition functions for each glueball state. At fixed volume, 
the pressure (= TlnZ/V) is then, 

P GB {T- e) = P^T- e) + P 2 (T; e) + ■ • • + P r (T; e) + • • • (27) 

The general term in Eq. (|2~7D corresponds to the pressure of a gas of 2D NR scalar r- 
particles at temperature T with masses M r given by Eq . (p3|) . That is, 

^^expj-^))^)" 2 ^ (28) 

An expected consequence of the previous result is that low T pressures vanish when we 
shrink the width strip to the line limit (recall M r (e) *— >°oo). However that is not the whole 
story. We can go beyond our qualitative expectation and obtain the small e behavior of the 
pressure exactly. If we make the e dependence in Eq.(|28|) explicit by using Eqs. (|23|) and 
(O), we find a non-analytic form for its e^O limit, 



InW^^^e-^ (-lne)^ ~? 7 = |(-) 1/2 

1 7T 



(29) 



5 Near the phase transition 

Pure gauge QCD on a narrow strip turns out to be a two phase theory. In section 3. we 
have studied its structure at very high T. We encountered a deconfining regime in which 
both transverse (0-particles) and longitudinal gluons (2D gluons) are the physical degrees 
of freedom. In contrast to the high T behavior, the Fock space of the theory at low T is 
made up of color singlet glueballs. Thermodynamics of the "cold" phase has no much to do 
to that of the free gluon plasma. Thermodynamic functions of the former correspond to an 
admixture of infinite species of glueball gases characterized by their different masses Eq.([27|). 
In the latter only transverse gluons, as a field gas, contribute to the partition function. 

We are going to investigate the critical behavior of the system both in its gluon plasma 
phase and in its glueball phase. If we are able to establish accurately the thermodynamic 
properties of the two phases in the neighborhood of the phase transition point, we will be in 
an optimal position to estimate its critical temperature. 

Let us begin analyzing the plasma phase. The starting point is the effective action for 
the gauge field modes Eq. (|i6|) , obtained after <p n integration. The partition function is then 
just, 



Z p iasma(T;e) = Z^(T;e) J[da n ]e s ef/ {a " ]/T 

Z^T-e) = J<1pEJ^Hp 2 +h 2 +^)] (go) 

Certainly the latter operation is not feasible in the most general case. However we 
will not need to work it out in an exact manner. Only very general properties of the 
functional integration will be required to understand the problem under consideration. The 
existence of two independent dimensionless parameters is the main ingredient to see why a 
complete knowledge of the effective action is not necessary. Because the theory possesses 
three parameters with dimension of mass (g2,n,T), we can form the dimensionless ratios 
(/i/T, g 2 / /i). Thus we can express any relevant function or functional in terms of them. For 



example, the high T calculation in section 3. is specified by the conditions (/i/T<Cl, g^//-**^ 
1). However nothing forces us to choose a priori fi<T, being both variables independent. 
We can take advantage of the freedom to select an arbitrarily large mass by working on a 
narrow enough strip. So that at any given T we can fulfill the conditions /i3>g 2 , /i3>T just 
by picking a small enough value of e. The first inequality is understood in the same way 
as in the calculation of the spectrum at zero temperature. That is the interaction term is 
small when compared to the kinetic piece of the bound state equation (NR approach). In 
this approach higher order contributions of the gauge field are suppressed as powers of g-ij \i- 
Therefore at any T we expect a weak coupling expansion of the effective action Eq.(|16|) in 
the gauge field modes to be valid when e<Cl. 

The structure of S e jf[a n ] in the weak coupling regime is easy to obtain from Eq.(|l6|) by 
expanding the logarithm in powers of (p 2 + fi 2 + u 2 )~ l {Mi + M2), 



S e ff[a n } = S der [d x a n ] + V e ff[a n ] 

V eff [a n ] = J dx\ l -ml{4) 2 + £ Xtn 2 nMAA 3 ) + - 
S der [d x a n ] = jdx^Y / (l + X' n )(d x a b n ) 2 + --^ (31) 

Although it is impossible to determine an explicit solution for all the coefficients, it is at 
our hand to establish their behavior for small e. We were able to calculate the first term of 
the effective potential in section 3. Although we only gave its high T expansion, it is also 
possible to show its form for n~^>T. The frequency sum in Eq.fll8|) can be transformed in 
an integral over the zero component of the momentum p || . The T-dependent part of this 
integral can be worked out, in the same way as the remaining p\ integration when we work in 
the /i>T regime. The whole operation results in a exponentially decaying behavior for the 
electric mass, ~e~ M /( 27rT ) e ^0. The exponential decay of the action coefficients is a general 
feature valid at every order for large enough masses. As a consequence all the coefficients of 
the different powers of a n appearing in Eq.([H]) vanish in the e— >0 limit. Moreover one finds 
that coefficients of higher powers decay more strongly than those of the lower ones. The 
above characteristic of the effective action coefficients is completely understandable within a 



3D framework. We know the 3D gauge theory tends to its 2D counterpart in the line limit. 
Thus all coefficients in the effective potential must disappear when e — leaving the pure 
QCD 2 action as the only remnant of color interaction. 

In summary. We find that at any given T, above the critical temperature, we can reduce 
the contribution from the effective action to that of a 2D gluon gas provided we take the 
strip width in a convenient way. Indeed the e — > limit is T-independent and therefore this 
limit is equivalent to write simultaneously (/i> g 2 , /i3>T) when T is taken at a fixed value. 

Mathematically we can express the above property as follows, according to Eqs. (|3~0l) and 

(0), 

Z plasma(T', e ) ~ -^(T; e) Z QC d 2 

Z QC d 2 = j[da n ]e-^£"[H(^K] (32) 

That is, the partition function of the gluon plasma factorizes in two pieces. One corre- 
sponds to a free gas of transverse gluons (Z^) and the other one to a free gas of longitudinal 
gluons (Zqcd 2 )- The QCD 2 gluon partition function is trivial. As mentioned in the intro- 
duction, it is T-independent and thus it has no thermodynamic significance as we will see in a 
moment. All the temperature dependence is produced by the quasi-free motion of transverse 
gluons. The plasma pressure is immediately obtained from Eq. (p2|) using standard finite T 
field theory techniques, 

PplasmaiT; e) ^ P^T; e) + P™ c = (iV c 2 - l)e~^ T ^ T 3/2 + P™ (33) 

Now comes the time to find out the properties of the glueball phase above the critical 
temperature. In section 3. we discovered the way the glueball admixture started to react 
under a small increase of the heat bath temperature. Only the lowest mass glueballs produced 
a discernible effect on the gas pressure, as Eq . (^7|) shows. The contributions to the global 
pressure of gases made out of heavy glueballs are exponentially suppressed at low T. In 
general, at any T the pressure of a given class of glueball is always exponentially suppressed 
respect to that of lighter boundstates. 



The transverse gluon is lighter than all its glueball boundstates. The pressure of the 
gas is always higher than any individual pressure of the glueball soup. Only the fact that we 
need to sum up a great amount of single contributions, from the infinite number of glueball 
states we can create from </> particles, can avoid that the total pressure of the glueball phase 
remains below the plasma pressure at any T. According to Gibbs criterion, if Pes < Pplasma 
at any T then the transition cannot occur. 

The mass distribution of the states becomes then crucial. The particular form of the 
density of the states determines the real shape of the pressure and thus the possibility for 
the phase transition to happen. One can really appreciate this point when writing the total 
pressure as 

P GB = jT 'dmp(m)e-^ T ^ T 3 / 2 (34) 

If the density of states increases rapidly for large masses (m ^> T) then it can compensate 
the thermodynamical suppression dictated by the Boltzmann factor and produce the neces- 
sary cut. The density of high energy states derived from the NR spectrum (B^) generates 
the following high T-dependence of the glueball pressure, 

p(m) m> ^ R m 1 / 2 P GB ~ T 7 / 2 > P plasma ~ T 3 / 2 at large T (35) 
As we can see, the transition occurs (see Fig. 2b). 

However although we can take for granted the existence of the deconfining transition, we 
are not describing properly the glueball phase near T c . The reason is we cannot extrapolate 
the NR spectrum (p^) to the highest energy boundstates. The accuracy of the calculation 
is based on the validity of the NR approach. For very energetical boundstates the binding 
energy is much bigger than p R . The typical NR parameter V^/ [j,r, an estimation of the 
characteristic relative velocity, is no longer small and a full relativistic approach has to be 
adopted instead. Moreover the heavy part of the spectrum has to contain necessarily multi-0 
states, not only the two particle sector. Pair excitations are now energetically available since 
typical momenta flowing through internal gluon lines are much bigger than the pair energy 
threshold. This means that our NR counting of glueball states is clearly underestimating 



the real number of them contributing to the total pressure. A careful counting of states is 
necessary in order to determine an accurate description of the glueball gas pressure near the 
critical temperature T c . 

Like in QCD 2 with fermions, one can use light cone quantization and obtain a full 
relativistic Bethe-Salpeter equation for multi-0 bound states valid for any constituent mass. 
This study has been already carried out quantitatively by Demeterfi et al. ||. In the 
description of the highest spectrum the constituent mass can be neglected respect to the 
boundstate mass [m ^> Thus we can study the growth of the density of the most 

massive states by analyzing their model in the limit in which the constituent mass equals 
zero. Their results in this limit are very interesting since they show that the density of states 



exhibits the Hagedorn behavior |10| , 



p H {m) « g b 2 - 1 am~ b e m/To [a] = 1 (36) 



Although only an estimation of the Hagedorn temperature T (e) = (1.4— 1.5) y N c /ir(g%/e) 
is given, the previous result is enough to obtain the relevant contribution to the glueball 
gas pressure near T c . The ultrarelativistic boundstates are much more abundant than the 
NR ones (compare the exponential behavior in Eq.(^) to the NR density). This results in 
a stronger raise of the glueball pressure with T and therefore in a decrease of the critical 
temperature (see Fig. 2a). The ultrarelativistic boundstates dominate the region below the 
phase transition point, as we will see next. 

In the pressure integral (p4|) it is legitimate to substitute the density of states by pu 
provided we work at energies above a typical mass M — A/i characterizing the ultrarelativistic 
regime. That is, A 3> 1. All states with masses bigger than M will be well described by 
Eq.©. 



P H (T;e) T *% e<<1 r dm p H {m)e~ m ' T T 3 / 2 + P™ 



~ (constant) ^( e )V^ e -^(^-D ( ZM _ + pvac 

(37) 

In principle, the constants Pl ac and Pqb m Eqs. (^) and (j37D may be different. The 



former corresponds to the pressure of the non-confining vacuum (free gluons), whereas the 
latter corresponds to the confining vacuum generating the Fock space of singlet color states. 
The difference between these two pressures is given by the difference between the vacuum 
thermodynamical potentials of the two phases (Q = —P) 



r 4> r GB — "non-conf "conf _ c conf 1 °° 1 

where £™£i * s the vacuum energy of the confining phase. 

The vacuum structure of pure gauge QCD 2 is trivial for it is a free massless theory, as it 
is apparent in a physical axial gauge (e.g. ai = 0). Obviously it is manifestly invariant under 
dilatation transformations and thus its vacuum energy is zero after zero point renormalization 
(£ vac ~<:T aa :>~<:d a jf:>=0). 

Pure gauge QCD% on a narrow strip is a confining theory. However the vacuum structure 



derived from the strip action (12) is the same as in pure gauge QCD 2 . After integrating the 
transverse fluctuation we obtain an effective action for the 2D gluon field whose first term 
is identical to that of pure gauge QCD 2 . This term renormalizes the coupling constant but 
in a momentum independent way, g 2 = Z(fi)g 2 . Thus no breaking of dilatation symmetry is 



encountered and we meet the same result as before [11 



Cnf ~<:^ ip :>~<:&(jf ) Strip :>= (39) 

The previous result can be likewise related to the absence of a gluon condensate (</ 2 > 
= 0) in the theory represented by the strip action (|I2"D . As we mentioned before the long 
distance effective action is proportional to that of pure gauge QCD 2 , which being a free 
massless theory cannot generate a non-zero gluon condensate. 

According to Gibbs criterion, at the critical temperature it must be true, (see Fig. 2) 



P GB (T C ; e) = P p i asma (T c ; e) (40) 

Once we know the two constants appearing in the expressions for the gluon and glueball 
gas pressure are the same, the equalization equation yielding T c is completely defined. The 



two constants cancel in Eq.(^0|) just rendering an equation for the T-dependent parts of Eqs. 



([33|) and fl37l) . We can now find the evolution of T c with the width valid for narrow strips. 

From the 3D point of view the two dimensional coupling constant is a width dependent 
object, (?2 = 9z/ e - If we wished to compare to numerical simulations in 3D, the natural 
variables to work with would be {g3,L) or (g 3 , e). The 3D coupling constant g% defines 
completely the color flux generated by the static quark- ant iquark pair in the infinite volume 
limit. The string tension, for example, can be expressed in terms of it since gf carries 
dimensions of mass. If we put the system on a strip (an asymmetric box, in the case of a 
lattice simulation), the natural question to answer is what happens to the infinite volume 
values of relevant observables such as the string tension, the glueball masses or the critical 
temperature under the new different boundary conditions. This is equivalent to fix the value 
of g% (maintain the infinite volume quark- ant iquark configuration) and tune only the strip 
width L. 

In Fig. 3 we give an approximate representation of the critical temperature as a function of 
e. It is obtained by solving numerically the equalization equation (|40|) for different values of 
e. We have observed that the discrete contribution Pqb is not relevant numericaly. However 
it is worth observing that the critical temperature does not exist in the exact line limit of the 
strip. This is not because the effective 2D coupling constant blows up in that limit, which 
pushes T c to infinity. In the strict e — > limit both the plasma and the glueball gas pressure 
are exactly zero independently of temperature. The equalization equation has no solution 
for e = 0. The previous feature is not related to the dimensional running of the 2D coupling 
constant, but to the effective decoupling of the transverse gluon /i# *— * oo in that limit. 
Therefore we recover in our strip approach a previously known result, pure gauge QCD 2 is 
always in the same (confining) phase. 

From the previous analysis one concludes that the existence of the two phase structure is 
a non-trivial quantum effect produced by the interaction between longitudinal and transverse 
degrees of freedom. Simpler models, as the non-compact 3D abelian gauge theory, lacks such 
property. As a consequence, they do not show the rich phase structure of the non-abelian 
theory p. 



6 Conclusions 



Our study establishes the preeminent role of transverse dynamics in the generation of a 
deconfining phase. Unlike their longitudinal counterparts, transverse gluons are dynamical 
and their dynamics is responsible for the formation of a screening mass for the 2D gluon. 
The loss of its massless character entails a radical change in a confinement order parameter 
such as the v.e.v. of the Polyakov loop. It becomes non-zero thus indicating the presence 
of a deconfining phase transition. It is highly remarkable that even a very small amount 
of transverse physics is enough to produce such a drastic change in the thermodynamics of 
pure gauge QCD 2 . 

Although the theory cannot be simple even for very small values of the strip width -it has 
to account for a transition phase-, an easier approach may be taken in that case. The Fock 
space of the theory changes radically depending on what side of the critical interphase we are 
describing. The partition function and all thermodynamic functions behave very differently 
below and above T c because of that circumstance. The advantage of the small e approach 
is it allows us to simplify considerably the dynamics of the physical states in both phases, 
but, at the same time, without spoiling the structure of their particular Fock spaces. In 
another words, in the small e regime we always find only color singlet states (glueballs) in 
the confining phase, and free gluon particles above T c . However both glueballs and gluons 
interact weakly among them in their respective phases. 

In this way we observe that the thermodynamics of the deconfining phase can be ac- 
curately described by a free gas of transverse and longitudinal gluons. Unquestionably we 
cannot predict the appearance of typical non-analyticities at the transition point with such 
a simple model. However the approach is as good as one wishes (just by manipulating e) 
provided we exclude the singular point. A clear signal that the non-analytical behavior is 
really included within the formalism is encountered in the effective action for the gauge field 
modes Eq . (|3TD . The non-zero modes are massless. If we integrated them in order to get the 
effective potential for the gauge field zero mode (an order parameter of confinement), we 
would face strong IR divergences already at one loop level. In finite temperature field theory 
this is a clear indication of a non-analytical behavior in T. This non-analyticity will show 



up in the form of infrared divergences in the perturbative expansion (even off mass-shell). 
Typical perturbative high T calculations break down, and one has to resort to complicated 
infinite resummations to find the non-analyticity behavior and predict T c fll2f . Fortunately 
in our case, with the exception of the critical point, all these effects are suppressed above T c 
when e is very small. 

The confining phase consists of weakly interacting gluon boundstates, which at low T 
behaves as a mix-up of inert glueball gasses. This fact is directly associated to the NR 
character of the low-lying spectrum. However the spectrum also contains fully relativistic 
states since there is no upper bound for the binding energy. These states form the high part 
of the spectrum. The high glueball spectrum is extremely populated because it is made out 
of many multi-gluon states whose energy levels are very close to each other. The number 
of such states grows fastly with the mass, as a relativistic treatment points out ||. They 
exhibit the Hagedorn behavior. 

After heating the system sufficiently we can get to excite the Hagedorn spectrum. How- 
ever this supply of heat does not turn immediately into an increasing of the gas temperature. 
Most of the energy is used to excite the great amount of massive states accessible to the 
system instead of increasing the thermal motion of the less energetical states. The energy 
necessary to produce a small change in the temperature becomes increasingly higher. The 
gas pressure blows up at a limiting temperature To(e) due to the exponential raise of the 
density of states. There exists then a temperature for which the pressure of the glueball gas 
becomes higher than that of the free gluon gas. Thermodynamically the latter is favored 
and the system evolves to this phase until the thermodynamical equilibrium is restored. 

From the microscopical point of view we can understand the phase transition as produced 
by a certain instability arising as a result of increasing fluctuations. We have seen that multi- 
gluon states are allowed by gauge dynamics on the strip. They are extremely heavy and thus 
the probability of a multi-gluon fluctuation is very small at low T. However we know that the 
pressure of the glueball gas grows rapidly near the Hagedorn temperature. The reason for 
such behavior is that an increasing amount of energy is needed to excite the large quantity 
of multi-gluon states with similar masses, formed by an increasing number of constituent 
gluons (this number grows with the glueball mass ||). Thus the probability of generating 



such multi-gluon fluctuations with a huge number of constituent gluons becomes bigger and 
bigger as we approach T . It is natural then to think of T c as the temperature for which 
the multi-gluon fluctuation involves an infinite number of gluons. The glueball gas turns 
unstable and the gluon plasma is formed. 

As it happened with the glueball spectrum of the model, the results here presented show 
an appealing agreement with our knowledge of the thermodynamics of higher dimensional 
gauge theories. It is interesting to remark why one expects these results to describe, at least 
qualitatively, the features of the 3D theory. Our approach consists in defining the whole 
3D theory under special boundary conditions (the strip compactification). The analysis 
developed here and in our previous work |jj accounts for the behavior of pure gauge QCD 3 
in the regime of very small witdths. The dynamics of the 2D effective action Eq . ([12]) , 
obtained to describe the long distance longitudinal physics of pure gauge QCD 3 , is defined 
completely by the 3D parameters (gs,L). The 2D effective coupling constant and the 2D 
scalar mass Eq.(|l3"l) are fixed by the 3D inputs. In this sense one deals with the real theory 
and not with a 2D model of it. In the latter we would take #2 and [i as free parameters to 
tune after comparing to 3D data. 

Certainly our analysis is restricted to a validity region (e<l) which lies far from the real- 
istic domain, the infinite width limit. However our unreachable non-perturbative transverse 
regime (e ~ 1) is the natural scenario for lattice simmulations. The study of asymmetric 
lattices (one lattice dimension substantially bigger than the others) complements perfectely 
our results. The small e regime is precisely the forbidden territory for the lattice simulation, 
which is limitated by the finite lattice spacing constraint. When both analysis merge, one 



discovers the two following relevant features [13 



1. One can extrapolate the small e results (i.e., in Eq. (|23|) ) to the e > 1 regime in a 
continuous way. This gives rise to a non-zero value for the QCD3 glueball masses and 
the string tension in the e^oo limit. 

2. There exists a stabilization of the e flow at a critical e* (corresponding to the critical 
width L*). The string tension and the glueball spectrum tend quickly to a constant, 
their QCD 3 value, above the critical point e*. 



The previous results strongly suggest that the description of the low energy longitudinal 
dynamics of pure gauge QCD 3 can be performed using a generalization of the 2D effective 
action Eq. flTjp . That is, through an effective action defined at the critical e*, 



^stripi^ ) Tr 



1 fS + (^) 2 + ^ + 9*2, e*) 



(41) 



.2' 

The low energy effective action contains only the terms with less number of covariant 
derivatives in the zero mode fields a a and <fi. All the parameters defining the 2D effective 
action get frozen at the critical width. In particular, the 2D coupling constant and the scalar 
mass acquire their infinite width value already at e*, (g2^>92-> yu. — //*) - The same happens 
to the effective potential U. 

The previous consideration allows us to understand why it is possible to get a sensible 
model of pure gauge QCD 3 using 2D adjoint matter gauge theories of the type (M) [|. 
In this approach, the coupling constant g<i and the scalar mass /x are no longer fixed by 
the 3D theory They are just input parameters. We cannot know them until we solve the 
underlying theory in the critical region. However it is reasonable to expect that there exist 
a set of 2D parameters which may describe accurately the low energy dynamics of the 3D 
theory. 

The qualitative features of the string tension and the glueball spectrum are kept along 
their trajectories in e space. They get stabilized at the critical value of the strip width L*, 
which is identified as the flux tube width. Therefore if the zero temperature properties of 
the model do not seem to disappear in the infinite volume limit, we have good reasons to 
believe that the qualitative picture here developed for the deconfinig phase transition can be 
also preserved. With maybe only one exception, the gluon condensate. Our analysis does 
not predict any difference between its values in both phases. 
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Figure Captions 



Fig.l: Diagrams contributing to the a n self-energy at high T. Dashed lines correspond to 
(f> modes, wavy lines correspond to gluon modes. 

Fig. 2 : Pressure versus temperature in the relativistic confining (dashed lines), NR con- 
fining (dot-dashed lines) and deconfining (solid lines) phase (g 3 — 1, e — 0.01). The vacuum 
pressure is taken to be zero in both phases (see text): (a) The critical crossing for a relativistic 
high energy spectrum; (b) The critical crossing for the NR spectrum. 

Fig3: Estimated evolution of the critical temperature with the strip width (L — e, g 3 = l). 



